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Abstract 

We present a local routing algorithm which guarantees delivery in all 
connected graphs embedded on a known surface of genus g. The algorithm 
transports 0{g\ogn) memory and finishes in time 0{g^n^), where n is the 
size of the graph. It requires access to a homology basis for the surface. 
This algorithm, GFR, may be viewed as a suitable generalization of Face 
Routing (FR) , the well-known algorithm for plane graphs, which we previ- 
ously showed does not guarantee delivery in graphs embedded on positive 
genus surfaces. The problem for such surfaces is the potential presence 
of homologically non-trivial closed walks which may be traversed by the 
right-hand rule. We use an interesting mathematical property of homol- 
ogy bases (proven in Lemma [4.3p to show that such walks will not impede 
GFR. FR is at the base of most routing algorithms used in modern (2D) 
ad hoc networks: these algorithms all involve additional local techniques 
to deal with edge-crossings so FR may be applied. GFR should be viewed 
in the same light, as a base algorithm which could for example be tailored 
to sensor networks on surfaces in 3D. Currently there are no known effi- 
cient local, logarithmic memory algorithms for 3D ad hoc networks. From 
a theoretical point of view our work suggests that the efficiency advan- 
tages from which FR benefits are related to the codimension one nature 
of an embedded graph in a surface rather than the flatness of that surface 
(planarity). 

1 Introduction 

Face routing (FR) was introduced by Kranakis, Urrutia et al in 1999 in 10 . It 
guarantees delivery in time linear in the size of the network, using logarithmic 
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memory, and was the first algorithm to exploit availability of position, now a 
common feature in many ad hoc networks. Its impact was twofold: it showed 
that it was possible to find local algorithms for ad hoc networks, as opposed to 
constantly re-adjusting routing tables, and it revealed the important role posi- 
tion could play in achieving this. The word "local" in this context means that 
the algorithm is executed by a mobile agent traveling through the network and 
using only local information at each node (including position) together with a 
"small" amount of transported memory (commonly restricted to be logarithmic 
in the size of the network) . 

More recently, in 2004, work of Reingold [M] showed universal exploration 
sequences (UXS's) are constructible in logspace thus implying that it is always 
theoretically possible to route locally in the above sense, even without position 
data (transporting memory of order O(logn) as overhead in the messages )0. 
However delivery time by this method is a high polynomial in the network size 
(not yet optimized), making it impractical for most ad hoc applications. This 
makes the critical challenge for routing in 3D ad hoc networks time efficiency 
and it opens the question of whether position on a surface can help improve this 
efficiency in 3D networks as it did in the 2D setting. We answer this question 
in the affirmative. 

Related recent work shows that, even using geometric information, no con- 
stant memory local routing algorithm exists for 3D UBW's that occupy a slab 
of thickness greater than (see Durocher et al |3j), however that work does 
not address algorithms using logarithmic memory. 

Another related result shows that no sublinear memory local routing algo- 
rithm exists for directed graphs, even if they are planar (Eraser et al [B]). This 
means no local routing algorithm exists at all in ad hoc networks in the pres- 
ence of uni-directional communication links, under the common (logarithmic- 
memory) definition of "local" used in the ad hoc community. 

In this paper we make a first step towards addressing local, logarithmic- 
memory routing algorithms for undirected 3D ad hoc networks, by considering 
those which lie on surfaces. In particular we are interested in ad hoc networks 
(for example of sensors) which may be deployed on terrains and outer surfaces 
of physical structures such as vehicles, buildings or other objects. We assume 
such surfaces are in fact oriented two-dimensional smooth manifolds which we 
will approximate by piecewise linear surfaces (PLS's). We discuss this in more 
detail in the next section. Given a PLS (and a homology basis for it) there exists 
a local, logarithmic-memory, quadratic-time routing algorithm, called general- 
ized Face Routing (GFR), which is correct for all connected graphs that are 
"reasonably" embedded on the surface. 

By "reasonable" we mean such that any two edges and/or reference curves 
(connecting curve, homology curves) may intersect each other only a bounded 
number of times (independent of graph size). We remark that Face Routing 
in the plane assumes edges of the graph to be straight lines and assumes a 
linear or piecewise linear connecting curve between sender and destination, thus 
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achieving this intersection property. Similarly in our case we take piecewise 
linear (PL) curves to guarantee reasonableness. 

Like Face Routing, GFR relies on the absence of edge-crossings. It is there- 
fore, like Face Routing, a base algorithm which would require further processing 
to handle edge-crossings and make it suitable for real- world ad hoc networks (see 
[UITT]). With appropriate assumptions on the network and surface it is possible 
to apply any of the standard methods for handling edge-crossings in the plane 
[l][TTJ[3] to surfaces, however the practicality of the assumptions involved would 
have to be studied in a more applied paper. We welcome input from readers 
with specific applications in mind. 

2 Terminology and Assumptions 

Assume Sg is a piecewise linear surface (PLS) of genus g that has been trian- 
gulated. This is a natural assumption in practice (see Section [2^ below) . We 
recall the genus of a surface is defined as the maximal number of disjoint, sim- 
ple closed curves which may be removed from the surface without disconnecting 
it. 

The algorithm GFR will be executed by an agent traveling in an embedded 
graph G on Eg. We assume, as is typical in ad hoc networks, that the time of 
transmission is small compared with the speed at which the network changes 
so that one may reasonably assume the network is fixed but unknown during 
execution. The agent does not have access to a global representation of the 
graph, which may be constantly changing, but only to the ID's of its current 
and neighboring vertices as well their positions on Hg. In addition it has access 
to a homology basis B for Eg . This will consist of closed simple curves contained 
in the 1-dimensional simplicial subcomplex of Eg . 

2.1 Homology and homology bases 

The reader who is not familiar with the concept of homology or homology ba- 
sis is referred to Jeff Erickson's [4J for an introduction to these concepts from 
a computational geometry perspective or to Hatcher [7j for a more detailed 
mathematical treatment. We give only a minimal overview here, of simpli- 
cial homology for surfaces. Assume a ring R. In this paper it will be Q. 
For k = 0,1,2, let Cfe(Eg;i?) be the set of all formal linear combinations of 
fc-simplices in Eg with coefficients in R. CkC^glR) is naturally a group under 
addition; its elements are called /c-chains. Define the boundary operator 

Sk ■ CkC^g', R) — > Cfe_i(Eg; i?) 

to be the linear map which sends each oriented simplex to the sum of of its ori- 
ented boundary simplices (or to zero if there are none). Those fc-chains which 
map to zero (i.e. have no boundary) are called k-cycles. They form a group 
which we denote Zk — ker6k for k = 1,2. For our case of i? = Q, connected 
1-cycles with integer coefficients can be interpreted as closed unparametrized 
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curves which are contained in the 1-dimensional simphcial sub-complex of 'Eg 
. Unconnected 1-cycles with integer coefficients are formal sums of connected 
ones. Let Bk = dk+i{Ck+i) denote the ^-boundaries, i.e. images of the bound- 
ary operator 5k+i for fc = 0, 1. Define the first homology group Hi{T,g] R) to be 
the quotient group Zi/Bi. Its elements are called homology classes. They 
are the equivalence classes for the equivalence relation "homologous to" : two 
1-cycles /3 and /3' are said to be homologous if their difference is a 1-boundary, 
i.e. /3 and — /?' are the two boundary components of some formal sum of 2- 
simplices. When i? = Q (or any other field), Hi{Y,g\ R) is in fact a vector space. 
A homology basis is technically a basis of this vector space and so consists of 
homology classes. But, in fact, one can always find a set of 1-cycles which are 
connected and have integer coefficients whose homology classes form a homology 
basis. As a result, it is common in the computational geometry literature to 
speak of the set of these closed curves (instead of their homology classes) as the 
homology basis. We will adopt this usage. We remark that iJi(Sg;Q) = Q^^ 
so a homology basis for Eg will consist of 2g closed curves. 

2.2 Closed curves 

As mentioned above, we may interpret 1-cycles with integer coefficients as closed 
(oriented) unparametrized curves or linear combinations thereof. Moreover all 
closed curves used in the homological arguments of this paper in Section [4] will 
be oriented and simplicial, i.e. are actually connected 1-cycles with integer 
coefficients. Therefore, in Section |4] we will write closed curve as a synonym 
for "connected 1-cycle with integer coefficients" , and will write simple closed 
curve if the curve has no repeated verticetH. Finally, when a parametrization 
is assumed, we will specify this and assume the orientation is compatible. See 
also Remark 13.41 

2.3 Planar Representation 

A set of disjoint simple closed curves fxi, . . . , fig on the surface T.g which do not 
disconnect Eg may always be completed to a homology basis 

The curves fii, . . . ^ fig also give a convenient way of representing the surface 
in the plane (see Fig. [T]lf|. Indeed, let C consist of the union of these curves. 

After removing C we know the resulting surface is homeomorphic to a multiply 

g 

punctured plane, namely some W = \ IJ Di, where the Di are disks. Once 

1=1 

we assume a particular homeomorphism $ from Eg \ C onto C , we can 

^The latter would unfortunately be a "cycle" in standard graph theoretic terminology; 
our convention avoids this. 

^For illustrations in the present paper we will additionally assume that /i^ n Xj is a single 
point ifi = j and empty otherwise (such a basis can always be found) , but this is not necessary 
to the algorithm; any homology basis will do. 
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map all points and curves of interest in Eg (that are not completely contained 
in the special C) to corresponding points and curves in the plane. We call this 
a planar representation of the surface. 

Moreover, an embedded graph G on T,g gives rise to ^{G), an embedded 
object (points and arcs) on U C K^, which we can interpret as a graph after 
making the identifications given by C: edges will be unions of curve segments 
with breaks occurring at elements of C. As the agent for GFR travels in G, all 
coordinates of neighbors it obtains will be given in terms of planar coordinates 
in U. See Section [231 for more details. 

We remark however that the proof of correctness that we give for GFR in 
Section 2] is not on the plane but rather on the PLS surface Sg , as this view is 
notationally simpler for the homological arguments involved. 



Figure 1: To the left: surface Eg of genus 5 = 3 viewed as the plane with 3 
handles attached (after removal of one point). To the right: Eg viewed as the 
plane with 6 holes pairwise connected by the A^. The set {fii, A^, : 1 < i < 3} is 
a homology basis for Eg. 



2.4 Note on Practical 3D Model and Future Work 

The property that Eg be a triangulated PLS arises naturally from the standard 
modeling schemes for representing real surfaces. We give two examples. 

Triangular Mesh Surface Models This is one of the most common model- 
ing schemes for representing an idealized surfaceQ in computer graphics. Such a 
presentation consists primarily of a set of triangles (given by their vertex coordi- 
nates in R^) and may also organize this information (e.g. in so-called "triangle 
strips") to facilitate stepping systematically from each triangle to a neighbor. 
Triangular Mesh Surfaces are in fact triangulated PLS's. Suppose we have fixed 
a three dimensional ambient coordinate system in which the real world object of 
interest (e.g. the surface of a building) is situated. Efficient algorithms exist for 
generating a PLS which approximates a smooth surface using samples from the 
surface (see for example [H]). One then triangulates the linear pieces if they 
are not already triangular so as to obtain a triangular mesh model. To take this 
a step further, suppose one has two triangular mesh surface models, Em and 
SouTj one on either side of the actual surface Eg and suitably close together. It 



triangular polytope) 
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is relatively straightforward to construct a piecewise linear map ^ from either 
Sin or EouT to the plane. If these two surfaces are built "compatibly" then 
one can uniquely map points lying between the two surfaces to points on the 
plane. This allows one to convert three dimensional coordinates of points on 
the real world Eg to coordinates of points on the plane. Also, straight lines 
connecting nearby points on the three dimensional surface Eg would be mapped 
to PL curves in the plane and one could then use standard planar techniques 
to deal with edge-crossings. The aspects of "compatibility" one would require 
of the triangle mesh models would depend on which edge-crossing techniques 
one wishes to use and also the particular setting. We welcome comments from 
people with particular sensor network settings in mind. 

Wire Frame Models Note that in structural analysis, wire frame models 
are often the primary objects used to model curved surfaces and an appropriate 
triangular mesh representation is then generated to fill out the wire frame model 
(see Om]). For simplicity, the union of homology curves Cm (resp. Cqut) could 
be taken to be contained in the original wire frame model. 

2.5 Data required by GFR 

We now summarize GFR's access to data. Each node has stored the model 
Eg (for example in form of a triangular mesh model) and the instructions for 
computing $ (see Section [2?3l) . The model Eg may be very approximate and the 
corresponding data small, depending on the application. Its main purpose is to 
record homological informatiord. It is moreover constant: it depends only on 
the surface and not on the number of nodes or their positions. Included in this 
data is the homology basis B. For example, assuming a triangular mesh model, 
this can be done by recording for every 1-simplex whether or not it belongs to 
the i'th curve of B (for i G ['2g]) and its position in a sequential numbering of 
the 1-simplices of that curve. 

Each node can (via GPS) obtain its momentary latitude, longitude and 
altitude. From these it computes (using $) its planar position on U C K^. 
The agent executing GFR can, at each node, query the neighboring nodes for 
their planar positions. 

When the agent is started by a sender S wishing to transmit to a destination 
T, it computes a connecting curve 7. This is any PL curve in the plane 
connecting $(S') to ^(r). 7 may be mostly composed of 1-simplices, to simplify 
calculations. 

GFR therefore operates essentially in the plane (within U C M."^) but uses 
identifications of certain curves in the plane and other homological data to 
accomplish its task. 



^ There exist various manifold reconstruction techniques from topological data analysis in 
recent years which yield good approximate manifolds with coinciding homological invariants. 
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3 Generalized Face Routing 

3.1 Basic Definitions and Face Routing 

Definition 3.1 (Border Walk). A walk P in G which has no edges on its left 
side will be called a border walk. 

Note that border walks are boundary components of regions, where a region 
is a connected component of the complement of G in Eg. In general, regions in Eg 
may have multiple boundary components. We assume the standard boundary 
orientation convention: a boundary curve (3 is oriented so that the region it 
bounds lies to the left when /3 is traversed in its preferred sense. Clearly each 
directed edge forms part of the boundary of exactly one region. We thus define: 

Definition 3.2 (Adjacent Border Walk). The (at most) two border walks which 
share a given (undirected) edge are said to be adjacent at that edge. 

Definition 3.3 (Trivial or Non Trivial Border Walks (NTBW), Tiled Region). 
// a border walk j3 is the only boundary component of some region, we say j3 
is a trivial border walk, otherwise it is a non-trivial border walk (abbreviated 
NTBW). In other words, a closed curve is an NTBW if and only if it is not 
homologous to zero. We define the tiled regior(^ Ti-{G) to be the union of G 
together with all regions of the graph which have only one oriented boundary 
component. 

remark 3.4. The graph G does not in general form part of the 1-dimensional 
subcomplex of Eg . Its vertices represent mobile nodes which may lie anywhere 
in Eg and do not necessarily coincide with vertices of the fixed triangulation we 
have assumed for Eg. G should for the time being be viewed as an embedded 
graph on Eg which gives rise to the object $(G) inU CM.'^. We do not specify 
how $ will do this. In Section \2.4\ we mentioned one way in which $ could be 
defined (on a neighborhood of Eg in M.^) so straight lines in between nodes 
would be mapped by $ to PL-curves in the plane. But many choices are possible 
depending on how one wishes to handle edge-crossings. Therefore we describe 
GFR in terms of the non-simplicial graph G embedded on Yig . 

In fact, for the proofs of Section^ we will actually subdivide the triangulation 
of Yig so that all curves considered are simplicial. But this is not done by the 
algorithm; it is only a convenience for our proofs. 

Regarding the (non-simplicial) definition of boundary orientation just given, 
we remark that it assumes the region in question has the same orientation as 
Tig. When we speak of 1-cycles bounding 2-chains, the latter may be formal sums 
with coefficients of any sign and 2-simplices themselves may have orientations 
that agree or disagree with that of Eg . 

In keeping with standard usage, we assume the right-hand rule assigns as 
exit edge at each node the next edge clockwise from the incoming edge. Such 

®We remark that we are departing from the definition of 7?.(G) in earUer papers by the 
author by here including in 7^(G) any region with a connected boundary, regardless of its 
genus. 
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a routing rule then traverses a border walk containing the starting (directed) 
edge. 

For each reference or connecting curve we assume an initial and final end- 
point and orient the curve from the former to the latter. For 7, these are S and 
T respectively, for a (closed) reference curve they both correspond to a single 
point on Y,g . Moreover, we assume a simple parametrization of each curve: to 
each point on the curve we associate its i-value, defined as the distance from 
the initial point to that point in the direction of orientation divided by the total 
length of the curve. This induces a total ordering of the crossings of any such 
curve with the edges of the graph (defining higher crossings as those with higher 
values) . 

We now define a slower variant of Face Routing. We assume a to be a 
connecting curve or homology curve and /3o to be a border walk that meets a 
(to start at S with a = 7, we let (3o be the border walk determined by the first 
exit edge clockwise from 7 at S*). In general, we start at a particular crossing 
a{to) of a with /Sq and for the purposes of the algorithm only, define next 
greater crossings in terms of t-values that have been shifted by subtracting 
^0 mod 1; note: we do not re-define actual values, just use shifted values to 
order crossings so as to make a{to) the least. 

Algorithm 3.5 (Modified Slow Face Routing (MSFR)). During the entire pro- 
cess, stop if T is encountered. 

1-13^ fio- 

2. Traverse j3 once to find a(tnext)> the next greater crossing with a after the 
entry crossing oi{tQ), keeping a record over all of (3 of the sum of crossings 
to the left minus crossings to the right for each homology curve. 

3. If any of the total sums is non-zero, stop and exit the algorithm, otherwise 
travel to a(tnoxt)- 

4-. Set j3 to the adjacent border walk, to <— tnoxt- Go to (2). 

Lemma 3.6. Let a be a connecting curve or a ^ B. MSFR fully follows each 
component of a Cl 7^(G) that it travels unless it meets T; more precisely, if 
an7^(G) starts at a{to) on the NTBW (3i and ends at a{t^) on the NTBW (32, 
then MSFR along a starting at a(io) on /3i stops at or T. And MSFR is 
reversible (performing MSFR on the oppositely oriented curve starting at a{t^) 
on /32 will surely reach j3i). MSFR along all components of a H 7^(G) takes a 
total time of 0{gn^) and memory 0{glogn), n= \G\. 

Although Face Routing finishes in linear time transporting logarithmic mem- 
ory, it does not have the first properties mentioned (following a within TZ{G) 
and reversibility). This is addressed in |5]. The time bound follows from the 
bound d on intersections of "reasonable" curves (since MSFR will travel the 
boundary of each region in TZ{G) as many times as that walk's edges meet a, 
the walk may have at most 0(n) edges and each edge may meet a at most d 
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times); at the same time, for each edge e traversed, MSFR must check all 2g 
curves of B for crossings with e. The factor of g in the space bound is likewise 
due to the total sums kept for all these curves (each total can reach be recorded 
with O(logn) bits). 

3.2 Generalized Face Routing (GFR) 

Let C consist of 7 together with all elements of B and all their oppositely oriented 
counterparts. Assume a fixed indexing of these 4^ + 1 elements beginning with 
7- 

For 7 and each curve a G B, we have already defined the t- value of a point 
on the curve; for the oppositely oriented counterparts oi a G B, we define the 
i- value analogously, reversing the roles of S and T. 

GFR transports a list of triples (z, ii, t9, where i is the index of a curve a in 
C, and ti,t^ are values along a. 

If MSFR along some a G £ stops before reaching T, then by Lemma [3.61 it 
must do so at a NTBW. If at this stage MSFR is begun on another a' G C, it 
will stop at another (possibly same) NTBW. We thus hop from one NTBW to 
another. 

Let r — (A/", £) be a virtual multi-graph whose vertices (A/") are all the 
NTBW's of the given embedded graph G, and whose edges (8) are connected 
curve- components c of a n TZ{G) : a G B U {7} where c C a is an edge 
between /3i and (32 , for /3i , /32 G A/" if and only if c has its endpoints at these 
NTBW's. An immediate consequence of this definition is: 

Lemma 3.7. Given an edge c = (/?i, /32) o/F, which corresponds to a component 
of a G BU{'y}, if MSFR along ±a is started on f3i at the corresponding endpoint 
of c then it will stop on (32- Hence, given any path /3i, /32, . • . , /3fe inT, such that 
T ^ S for any border walk 5 meeting one of the ct = {/3i,/3i+i) : i — 1, . . . ,k — 1, 
reverse MSFR can be used iteratively to return to /3i from j3k as long as there 
is available a list of the t-values corresponding to the initial and final endpoints 
of the Ci, together with the index of the curve in C to which each Ci belongs. 

Using this lemma, we define: 

Algorithm 3.8 (Generalized Face Routing (GFR)). During the entire process, 
exit if T is encountered. 

1. Use MSFR on 7 starting from S until it stops. The agent is now at a 
NTBW, i.e. vertex o/F. 

2. Perform a depth first search o/F, keeping - at each subsequent NTBW /? 
where MSFR stops - a record of {i,ti,t'j), where i is the index of the curve 

^If a has no self-intersections then a connected curve-like component is just a connected 
component. If q does have self-intersections, we assume it is parametrized so q is a home- 
omorphic image of the unit circle. Then retain only that part V of the circle mapped into 
'1Z(G). Connected curve-like components of ar\'JZ{G) are images of the connected components 
of P. 
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ai € C just followed, ti is the t-value along a where MSFR was initiated 
and t[ is the t-value along a where MSFR reached /3. 

remark 3.9. That an NTBW has been previously visited is detected by GFR 
by the presence of a crossing ai{ti) such that {i^ti^t[) and ai correspond, as 
described, to one of the recorded triples. 

In the next Section, we will prove: 

Proposition 3.10. There exists a path in T from first^j) to last{'j), resp. the 
first and last NTBW's met by 7. Moreover, \T\ < 2g. 

This implies GFR must eventually reach last{'j), from which MSFR will be 
performed on 7 until T is reached, thus proving the following Theorem. Fig. [5] 
illustrates the process. 




Figure 2: This is a planar representation of a surface Eg (with g = 4 in this 
example). The agent executing GFR sees local portions only of this picture. 
The shaded area is (the image under $ of) 7^(G) - we will omit writing $ and 
keep all notation from E. The curves /3i,/32,/33 are three of the NTBW's (two 
more are unlabeled). The dotted piecewise linear curve is 7; i.e., the connecting 
curve between S and T. The pink disks are the closed disks bounded by 
and —fli respectively; the thin solid curves connecting pink discs in pairs are 
the Ai {fii paired with —fLi). GFR on these data follows the arrows in boldface, 
recording a triple at each star in the numbered sequence. Note that appears 
in two pieces in the plane, but is in reality a single curve (when pink disk pairs 
are identified). Also, in going from the second star to the third, GFR will first 
take the arrow into the annular shaded area and then re-emerge in the central 
disk-like shaded area before arriving at /^s . From /Js the curve 7 is used to reach 
T. 
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Theorem 3.11. GFR is a local routing algorithm transporting memory of 
0{g\ogn), and guaranteeing delivery in time 0{g^n^). 

The time bound follows since |r| < 2g and also there are 0{g) reference 
curves in C while all edges of F which correspond to any one of these curves 
can be traversed by MSFR in time 0{gn^) by Proposition 13.101 The record 
of triples kept by GFR has size 0{glogn) and traversing each edge of F also 
requires memory 0{g log n) by the Lemma, hence the space bound. 

4 Proof of Proposition 

We now prove Proposition 13.101 It relies on a mathematical fact, Lemma [4.3( 
which we prove separately below. From now on we assume that the original 
triangulation of Eg has been subdivided so that all closed curves referred to 
- including those which are NTBW's - are simplicial. We do this to prove 
Proposition only. The algorithm does not do this in practice. 

Proof. We divide the proof into two parts; showing first^-y) and last{'-f) are in 
the same connected component of F, and establishing the bound on the order 
of F. 

(1) As we follow 7 from first{j) to last{j), let first{-f) — /32, ■ • • , /3n-i, /?« — 
last(j) be the NTBW's encountered. Successive NTBW's Pi^Pi+i are alternat- 
ingly either: 

- connected to each other by a component of 7 n 7?.(G) or are 

- boundary components of a common region F . 

In the first case, there is by definition an edge between /3i and Pi+i in F. In the 
second case, by Lemma 14.31 there is a path in F between them. Thus we obtain 
a path in F from first{'y) to last^j). 

(2) To count the NTBW's, suppose there are N of them and let T be the 
collection of all "non-trivial" regions, meaning regions which have more than 
one boundary component. Let k be the number of elements in J-'. We may cut 
along all but one NTBW per element of and obtain a surface E which is still 
connected, since each region that is an element of will still be attached to 
71(0), which itself is connected. We are cutting along N—k curves, so N — k < g. 
On the other hand, cutting along just one NTBW of each element of T would 
not disconnect the surface either, so k < g. Hence, N < g + k < 2g. □ 

remark 4.1. Recall that part (1) of the proof of Proposition \3.lO\ was divided 
into two cases. In the first, we considered vertices of T (i.e. NTBW's of G) 
which were adjacent in F since the NTBW's involved were connected by pieces 
0/ 7 n TZ{G). In the second case we appealed to Lemma \4.3\ to establish a path 
in F between two NTBW's which share a common region F. In fact such a 
path exists without using any pieces of y D Ti-iG) . In other words if F* is the 
subgraph ofT obtained by removing all edges of the form ^OTZ{G) then the path 
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of interest exists even in T* . We state and prove the Lemma in this stronger 
form. 

remark 4.2. In the Lemma we will also allow B to contain formal sums of 
simple closed curves, since the added generality simplifies the proof by induction. 

Lemma 4.3. [f| Given a set B of 1-cycles with integer coefficients generating 
iJi(Eg;Q) and a connected graph G embedded in Y.g, define T* to be a virtual 
multi-graph, with the same vertices as T in Section \3.S\ but with edges being 
connected curve- components of aC\TZ{G) : a ^ B only. Then whenever two 
NTBW's 13 and /?' are distinct boundary components of the same region ("face") 
F determined by G, there exists a path in T* from j3 to 13' . 

4.1 Intersection Number 

In order to prove Lemma 14.31 we will use certain facts about homology bases 
expressed in terms of the notion of intersection number. We therefore define 
intersection number and list its relevant properties, sketching their proof in 
elementary terms. A formal treatment may be found in [8 . We recall that is a 
simplicial complex and we use the term closed curve to mean a connected 1-cycle 
(we are using simplicial homology; see the brief introduction in Section |2.1[) . 
Every closed curve may be written as a sum of simple closed curves, i.e. 1- 
cycles with no self-intersections. Now, given any two simple closed curves a, 
we define their intersection number ^{a,l3). If a, 13 coincide, we set /3) 
to zero. If they do not coincide, let tr be a shared connected 1-chain that is 
maximal in the sense that at its two distinct endpoints, qi and (j2, ct and /? 
diverge. We assume that qi is the start of the shared chain and q2 the end, 
using the orientation of a (which may not coincide with that of j3). On the 
other hand, we use the orientation of /3 to establish a left and a right side of fi 
(here the orientation of Eg is implicit as well). Now define L{a) to be zero if a is 
incident to the same side of f3 at qi and (72, and define L{a) to be +1 (resp. —1) 
if a is to the left of (3 at qi and then to the right at 52 (resp. right and then left). 
Positive t(cr) amounts to saying a leaves a as an outward normal to the region 
bounded by (3. Now, suppose this has been done for all shared sequences a and 
sum i{a) over them to obtain ^{a,P). This defines intersection number for all 
simple closed curves. Extend this definition by bilinearity to formal linear sums 
of simple closed curves. The key properties of intersection number that we will 
use are: 

1. #(a, (3) is a homology invariant, i.e. depends only on the homology classes 

*This may be stated more generally, without reference to a graph: given TZ, a connected 
subsurface with boundary of a Riemann surface Sg with fixed choice of homology basis B, 
and two closed curves /3,/3' which are boundary components of a single component of TZ'^, 
there exists a sequence of homologically non-trivial curves /3 = /3i , /32 , ■ • • , 1 , /3n = /3' such 
that each pair of successive curves I3i and is connected by a curve segment in the family 

{antZ : aeB}. 
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2. ^{a, [3) is bilinear 

3. #(a,/3) is non-degenerate, i.e. any curve which is not homologous to 
zero must have nonzero intersection number with one of the elements of a 
homology basis. 

These properties show that intersection number provides a non degenerate bilin- 
ear form on the first homology vector space, and they express Poincare duality 
in the special case of oriented surfaces (see [TT] for an accessible treatment of 
Poincare duality and intersection forms). We outline a proof of the listed prop- 
erties alone, giving the essential ideas. 

Proof sketch, [of properties above] We now assume all curves are simple (i.e. 
have no self-intersections). The general case follows by straightforward exten- 
sion (property 2. immediately), since arbitrary closed curves are formal linear 
combinations of simple ones. Moreover, for ease of discussion, we parametrize 
each curve compatibly with its orientation so we have an order in which to follow 
the curve through all its 1-simplices. 

For property 1., note that if (3 and /3 are homologous then /3 and — /3 are the 
two boundary components of some region F. As mentioned in the definition, 
after a shared chain a with positive L{a), a exits the region F, while after a 
shared chain with negative t(p), a enters F. Since any entry (resp. exit) across 
/3 that does not have a corresponding exit (resp. entry) across (3 must have 
one across (3, we conclude that ^{a, (3) = jj={a, (3). The argument in the first 
coordinate is analogous. Property 2. follows by a similar basic argument since, 
for example, having (3 — (3' + (3" implies — /? and /?' and /3" are the boundary 
components of a region F. Finally for property 3., note that if ^{a, f3) = for 
all a in a homology basis B, then by property 2. /3 would have zero intersection 
with all curves of the surface. If f3 were not homologous to zero then removing 
/3 from Eg would not disconnect and so there would exist a path within the 
1-skeleton from an edge on one side of (3 to an edge on the other side of (3, not 
crossing P itself. But now, by connecting the two endpoints of this path (to (3 
and possibly along it) one would obtain a curve which has intersection number 
±1 with /?, a contradiction. □ 

Definition 4.4 (Crossing). We will refer to a maximal shared connected 1-chain 
of a and 13 as a crossing whenever a and j3 are on opposite sides of each other 
at the start and end of the chain. We will say a crosses /3 if such a crossing 
exists (even if the intersection number of the two curves is zero). 

remark 4.5. Any 1-cycle rj which has zero intersection number with a closed 
curve (3 is homologous to a sum of 1- cycles which can be parametrized (compat- 
ible with orientation) so they do not cross /3. Indeed, one pairs crossings in one 
direction with crossings in the other and then connects a pair of crossings by 
two copies of a portion of (3, with suitable (opposite) orientations; when taking 
the formal sum, these portions of (3 cancel out. 

On the other hand, for any 1-cycle rj that has non- zero intersection number 
with a closed curve f3 there must exist a curve-lik^ connected component c ofri\f3 
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which starts on one side of /? and ends on the other. Indeed if this were not so 
then following rj (according to any fixed orientation- compatible parametrization) 
the values of i at the crossings with 13 would alternate between ±1, thus giving 
4f^{ri,j3) —0, a contradiction. 

4.2 Proof of Mathematical Lemma 

of Lemma We prove the Lemma by induction on the number of vertices 
of r*. For or 1 vertices it is trivially true. Now suppose F* has exactly two 
vertices, which correspond to the NTBW's /3 and /?'. Then these two NTBW's 
must bound the same region F (since an NTWB by definition is a boundary 
component of a region with multiple boundary components). By property 3. of 
intersection number there exists a € B such that ^{a, /3) ^ 0. Therefore by 
Remark l4 . 5 1 there exists at least one component of There is moreover at least one 
such component c. So /3 and /?' are connected by c and hence F* is connected. 

Now assume the statement of the Lemma holds whenever F* has k or fewer 
vertices. Suppose F* has fc + 1 vertices and k > 2. Let (3 and /3' be two distinct 
NTBW's which bound a common region F, and let /3" be a third NTBW, 
distinct from /3 and /?'. Now suppose we "cap off' (3". More specifically, by 
capping off we mean cutting along f3" and gluing in a small triangulated PLS 
surface homcomorphic to a closed disk on either side of the cut. The resulting 
surface E' has genus g — 1 and it remains connected since G is (moreover, it is 
PLS and triangulated). The caps we glue in to /3" and —/3" are both closed disks 
and we will retain the name /3" for the boundary of the disk capping off /3" . 
Moreover, since f3" was a subgraph of G - in fact the only part of G that was 
affected by the cutting and gluing - let the new f3" still form the same part of G 
and continue to use the name G for the resulting embedded graph on S'. Notice 
that for this graph the tiled region will be extended by exactly one disk and f3" 
will no longer be a NTBW but all other NTBW's in Sg wiU stiU be NTBW's 
in E'. This is because even if /3" in Eg was a boundary component of F, there 
would still remain the two boundary components (3 and /?' for the region after 
capping off and by definition of NTBW, both /3 and l3' would still be NTBW's. 
If we denote the new tiled region in E' by 7^(G)', we have 7^(G) C TZ{Gy . 

By a standard homological argument (given in Lemma [4.6l below) there exists 
a a set B' of 1-cycles with integer coefficients generating _ffi(E';(Q)) such that 
for each a G B', every connected curve-like component of a fl Tl{G)' is either 
a connected curve-like component of ?7 fl 'R-{G) for some 77 G S, or else a union 
of portions of /3" together with connected curve-like components of 77 n 7?,(G) 
that meet f3" . If we define (F*)' to be the corresponding virtual multigraph, 
with vertices that are NTBW's in E' and edges that are connected curve-like 
components of a n 7?.(G)' for a ^ B' , then we see that the vertex set of (F*)' 
is obtained from that of F* by eliminating the single vertex f3" (and no other). 
Moreover (using this correspondence of vertices) the property of B' given above 
means that any edge of (F*)' (from v to v' say) is either an edge of F* or one 
can find two edges in F*: one from v to /3", the other from f3" to v' . Thus a 
path between vertices in (F*)' implies a path between those same vertices in F*. 
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We are now done since we know by the inductive hypothesis that there is a 
path connecting /? and /?' in (F*)'. □ 





Figure 3: Induction on the genus. The top surface is Sg, the bottom E' obtained 
by cutting along the NTBW first{"i) and gluing a disk in on each side. The 
darkly shaded area represents the tiled region in each surface. The red and 
green curves in the top surface are its \i and Hi respectively (same surface as 
in Figure [J). The red curve in the middle plays the role of A, by means of 
which we define the new reference curves for S' as shown on the bottom surface 
(red/green retained). If we solve the routing problem based on a connecting 
curve 9 in the surface S' (which has been taken to be A \ first{'y) and is shown 
dashed) this will give us a solution to the original problem which was in terms 
of 7 in Sg. But the surface E' has genus strictly less than g, so we are done by 
induction. 

The argument we just made is illustrated in Fig. [3] in terms of routing. We 
now state and prove the result we used concerning homology bases and capping 
off. 

Lemma 4.6. Let B be a set of l-cycles with integer coefficients which generates 
HiCEg-jQ) and T,' the surface obtained by capping off a closed homologically 
non-trivial curve /3 (and — /3 too). Then there exists another set B' of l-cycles 
with integer coefficients which generates iJi(E';Q) such that for each a £ B' , 
every connected curve- lik^ component of ar\TZ{G)' is either a connected curve- 
like component of rj Cl 'R-{G) for some t] £ B, or else a union of portions of j3 
together with connected curve-like components of rj r\TZ{G) that meet 13. 

Proof. We construct a set of l-cycles which generates the homology of E' and 
has the required properties. 
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Since /3 is not homologous to zero, there is some a € B such that #(a, /?) 7^ 
by property 3. above. Denote it A and let m = #(A, (3) ^ 0. Now consider the 
following set of 1-cycles with integer coefficients: 

5 = {ma-#(a,^)A : aeB}. 

All of these curves have zero intersection number with f3. Using Remark I4.5( 
each one is homologous to a sum of (parametrized) closed curves not crossing 
(3. Let S' be the set of all these closed curves; they are also closed curves in E'. 
We claim the elements of S' generate the first homology of S'. 

Warning: when we write [•] we mean a homology class in Hi(T,g,Q); the 
homology classes in Hi{Y,g,Q) we will simply describe as such and not give 
them a short-hand notation. 

Let S be an arbitrary closed curve in S'. Since /3 bounds an open region D 
in E' that is homeomorphic to an open disk, any piece of S which is a chord of 
D is homologous (in S') to a connected sub-chain of (3. Using this to substitute 
all such chords of S, we obtain that S is homologous (in E') to a curve which 
does not cross (3 and in fact lies outside D. Without loss of generality we may 
take this new curve as 6 (since we are interested in generating the homology of 
E'). It corresponds to a curve in Eg (by replacing 1-simplices on the boundary 
of the cap with 1-simplices of /3, as usual). Denote this curve in Eg also by S. 
It is homologous in Eg to a linear combination of the elements of B, say 



with coefficients a{a) £ Q indexed hy a G B. We have 
= #(<5,/3) = ^a(a)#(a,/3). 



Therefore, 

' ^ m 



^ ^ ^ [ma] «(«)#(«, /3) [A] 



None of the elements of S' crosses /3 and so neither does a linear combination 
of them. We thus have two curves, ^ {ma — ^{a,f3)\} and S which lie 

aeB ™ 

outside D and are homologous to each other in Eg. Denote the first curve 77 
and let k be a 2-chain in Eg with boundary 77 — (5. When we cap off /? as 
described earlier we produce the curve ry — 5 in E and a new 2-chain k' in E^ 
whose boundary is rj — S. Thus rj and S are also homologous in E'. So S' indeed 
generates the first homology of E'. □ 
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5 Conclusions 

Face Routing (FR) is at the base of most position-based routing algorithms used 
in ad hoc networks today. It is a logarithmic-memory, local algorithm guarantee- 
ing delivery in embedded graphs on the plane but not on positive genus surfaces. 
We have exhibited a position-based algorithm, GFR, which guarantees delivery 
for embedded graphs on surfaces of arbitrary genus. It is also local and uses 
logarithmic memory. Like FR, it is a base algorithm which as a next step could 
be adapted to handle edge-crossings. We welcome communication on particular 
settings of interest, as the method of handling edge-crossings would have to be 
somewhat tailored to the application. While universal exploration sequences 
(UXS's) provide a non position-based logarithmic-memory, local routing algo- 
rithm for any network, the time is a high polynomial (> 16) and no better 
algorithm has until now been proposed which would apply to 3D ad hoc net- 
works. FR and GFR take linear and quadratic time respectively. By restricting 
to the setting of graphs on surfaces which retains the codimension-one character 
of graphs in the plane, we have obtained (with GFR) an FR-like gain in time 
efficiency for the task of local, logarithmic-memory routing in a class of 3D ad 
hoc network. 
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